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ABSTRACT 

We report the detection of a massive rotating disk around the high-mass 
Class candidate NGC 7538 S. The disk is well-resolved with BIMA (9 A =3"7) 
in 3.4 mm continuum and in H 13 CN J=l— >0. It is seen nearly edge on and has 
a size of ~30,000 AU. A young, powerful outflow perpendicular to the rotating 
disk is mapped in SiO J=2— >1 and HC0 + J=l— >0. The dynamical age of the 
outflow is < 10,000 yr. The velocity gradient seen in H 13 CN is consistent with 
Keplerian rotation. Assuming that the gas is gravitationally bound, the mass of 
the central object is ~ 40 M Q . The mass of the continuum "disk" is > 100 M 
and has a luminosity of ~10 4 L . H 13 CN gives a mass ~400 M for the rotating 
disk, and ~1000 M for the extended (20") envelope. Our observations confirm 
that this is an extremely massive protostar in its earliest stages. 

Subject headings: ISM: clouds - (stars:) circumstellar matter - stars: formation 
- stars: pre-main sequence - submillimeter 

1. Introduction 

The formation of high mass stars is still poorly understood, although one expects a 
high mass protostar to form an accretion disk and drive an outflow similar to what one sees 
in low mass protostars. Recent surveys show that outflows are very common in high-mass 
star forming regions (Beuther et al. 2002, and references therein) but that accretion disks 
have been far more elusive, even though they will have to be larger and have larger velocity 
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gradients than accretion disks around low mass stars, which show velocity gradients of only 
a few tenths of km s" 1 , see e.g. Mannings and Sargent (1997). Although there have been 
a number of papers reporting rotating disks around high-mass stars, most of these show 
marginal, barely resolved systems. The best case to date, the young high-mass star IRAS 
20126+4104 (Zhang et al. 1998; Cesaroni et al. 1997), has a bolometric luminosity of 1.3 10 4 
L and drives a massive molecular outflow. 

In this letter we present the first results from high spatial resolution BIMA observations 
of NGC7538S. Sandell and Sievers (2003) showed that NGC7538S is a cold, very massive 
elliptical dust source (~4OOM ) associated with H2O and OH maser emission but only weak 
free-free emission, which has all the characteristics of a high-mass Class source. The proto- 
stellar source is located in the molecular cloud south-east of the large H II region NGC 7538, 
which is known to contain several centers of active and on-going high-mass star formation 
(Werner et al. 1979). The distance to NGC 7538 is assumed to be 2.8 kpc. NGC 7538 S is 
about 80" to the south of the well studied ultracompact H II region IRS 1. 

NGC 7538 S is more isolated than IRS 1 and therefore easier to study. We confirm 
that NGC 7538 S is a high mass protostar and resolve the disk around the protostar both 
in continuum and in molecular lines. In this letter we present data that suggest that the 
protostar is surrounded by a massive, rotating accretion disk, which drives a very young and 
powerful outflow perpendicular to the disk. 



2. Observations and Data Reduction 

The observations of NGC 7538 S were made with the BIMA array in 2001/2002 using 
two frequency settings in the B and C-array configuration. The correlator was split into 
four 25 MHz bands resulting in a velocity resolution of ~0.34 km s -1 . In this letter we base 
most of our discussion on the frequency setting that included H 13 CN J = l— >-0, and HCO + 
J=l— »0, which was observed in both B- and C- arrays. The other frequency setting included 
the molecular transitions HN 13 C J=l->0, SiO v=0 J = 2^1 and H 13 CO+ J=l->0. Only 
one C-array track was obtained in this setting. The images therefore have poorer spatial 
resolution and sensitivity, but are consistent with the deeper images seen in H 13 CN and 
HCO+. 

The data were reduced and imaged in a standard way using MIRIAD software (Sault, 
Teuben and Wright 1995). Phase calibration was applied using observations of the quasar 
0102+584 at intervals of 30 minutes. The phase calibrator was observed using an 800 MHz 
bandwidth for 3 min. 3C454.3 was observed for 10 minutes as a bandpass calibration. The 
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flux density scale was checked from observations of Mars. The uncertainty in the absolute 
amplitude scale ~15 %, but the relative amplitude of spectral lines within the same receiver 
tuning is within ~5 %. 

The data were imaged with weighting inversely proportional to the variance in order to 
obtain the best signal to noise ratio. Averaging 4 spectral channels we obtained an RMS noise 
level of 24 mJy (0.3 K) with a synthesized beam FWHM 3"75 x 3"59 and peak sidelobe level 
of -7 %. Spectral windows which did not contain any significant spectral line emission were 
averaged to provide a continuum image with an RMS ~3 mJy. The continuum emission was 
subtracted from the spectral line channels, and the images deconvolved using the CLEAN 
algorithm. 

3. The continuum "disk" 

The single dish sub-mm continuum maps (Sandell and Sievers 2003) show NGC 7538 S 
as an elliptical source embedded in a narrow dust filament. With BIMA we filter out the 
extended cloud emission and detect only the sub-mm source. At 3.4 mm the continuum 
source is well resolved with a peak on the northeastern side of an extended elliptical source 
(Fig 1). A two-component Gaussian fit shows that the northeastern peak is an unresolved 
point-like source with a flux density of 20 mJy offset 3"7, 0"5 from the nominal sub-mm 
position. The position of the 3.4 mm point source coincides with the OH and H 2 maser 
and the 6 cm VLA position (Argon, Reid and Menten 2000; Kameya et al. 1990, O. Kameya 
2003, private communication). At 6 cm the VLA source has a flux density of 2.8 mJy. If the 
6 cm emission is due to an ionized wind, i.e. follows a z/ 6 frequency dependence, about half 
of the 3.4 mm flux is an excess, presumably due to hot dust close to the protostar. However, 
the free-free emission could also have a steeper frequency dependence. In the following we 
assume that the point-source emission is all due to free-free emission. 

The extended elliptical source has a size of 13''7 x 8''1 and a position angle (P.A) of 71°, 
offset +V/7, —O'.'A from the nominal sub-mm position with an integrated flux density of 190 
mJy, i.e. it agrees within errors both in size and position with the sub-millimeter source. In 
this respect this source differs from IRS 1 and IRS 9, which both appear unresolved at 3 mm 
(van der Tak et al. 2000) and from nearby low mass protostars, which are only marginally 
resolved with a ~0"5 beam (Looney, Mundy and Welch 2000). 

To estimate a mass from our continuum observations we adopt a single-temperature 
optically thin thermal dust model. The total mass of gas and dust, M, can then be expressed 
as M = Sj,D 2 /(Kj,Lv(T d )), where B„(T d ) is the Planck function, k v is the dust mass opacity, 
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Td is the dust temperature, D is the distance, and S u is the integrated flux density at the 
frequency v. We assume a gas to dust ratio of 100 and adopt a dust mass opacity at 87 GHz, 
k 87 = 0.0073 cm 2 g~ 1 , corresponding to a dust emissivity, (3—1. This value was also used 
by Looney, Mundy and Welch (2000) in their BIMA study of low mass protostellar disks. 
With a dust temperature, = 35 K Sandell and Sievers (2003), we find the total mass of 
the continuum disk to be 100 M . We note, however, that the isothermal greybody fits by 
Sandell and Sievers (2003) predicted a steeper emissivity law, (3 = 1.6, which would result 
in a five times smaller mass opacity and hence a more massive disk. 



4. A rotating protostellar disk ? 

The map of integrated H 13 CN J=1^0 emission (Fig. 1) shows a bright elongated 
source surrounded by a more extended envelope superposed on the 3.4 mm continuum. The 
elongated source in the center of the map has the same orientation and extent as the source 
we see in dust continuum. For a two-component Gaussian fit, we find the bright emission 
offset 2''0, 1"0 from the nominal sub-mm position, i.e. within errors coincident with the 
6 cm VLA and the OH/H 2 maser position at +3"7, + 1"0. The fitted size: 11"1 x 6"8 
with a P.A. = 55 ± 5°, is slightly smaller than the continuum disk. The position velocity 
plot taken along the major axis of the bright H 13 CN emission shows a clear velocity gradient 
across the source (Fig 2). The emission is red-shifted to the southwest and blue-shifted to 
the northeast with the center of symmetry coinciding within 2" with the position of the 
protostar. A rotation curve derived from these data is slightly asymmetric with the central 
velocity close to —56 km s -1 . It reaches a peak of ±1.35 km s _1 at a radius of 5" (14000 au) 
in both directions. The observed velocity gradient corresponds to an enclosed mass of ~ 30 
/ sin 2 i M Q for a rotating disk to be gravitationally bound. With an assumed inclination, i 
= 60° (see below) the mass is ~ 40 M . 

H 13 CO + J = 1 — ^ and HN 13 C J = 1 — ^ were observed with only one C-array track 
and therefore lack the spatial resolution to resolve the disk. Both molecules show strong 
extended emission towards the protostar and we can therefore use all three molecules to get 
an estimate of the mass of the disk and the surrounding molecular envelope. For estimating 
column densities and masses we assume LTE and that the excitation temperature for all 
molecules is the same as that for the dust, i.e. 35 K. The permanent dipole moments are 
taken from Blake et al. (1987), except for HCO + , for which we adopt fi = 4.07 Debye (Haese 
and Woods 1979). For total mass estimates we use abundance ratios to H2 similar to those 
in the OMC-1 extended ridge (Blake et al. 1987), i.e. [HCO+]/[H 2 ] = 2 10~ 9 , [HCN]/[H 2 ] 
= 5 10- 9 , and [HNC]/[H 2 ] = 5 HT 10 , and the isotope ratio, [ 12 C]/[ 13 C] = 85 (Wilson and 
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Rood 1994). For H 13 CN we get a peak column density of 7 10 13 cm~ 2 for the disk and 5 10 13 
cm -2 for the surrounding envelope, resulting in a disk mass ~ 400 M and an envelope ~ 
1000 M Q . For H 13 CO + and HN 13 C we do not have enough spatial resolution to separate the 
disk and the envelope, but we derive combined masses ~ 1000 and 3000 M Q , for H 13 CO + 
and HN 13 C, respectively. 



5. The high velocity outflow 

We see clear high velocity wings in both HCO + J=l — ^0 and SiO v=0 J=2— >1. Towards 
the protostar the HCO + shows dominantly blue-shifted high velocity emission extending 29 
km s _1 from the systemic velocity (~ —56 km with a sharp, deep absorption on the red- 
side at nearby cloud velocities, suggesting infall of the surrounding cooler envelope. The SiO 
line wings at the position of the disk appear more symmetric, presumably due to the broader 
beam and the absence of self-absorption. Both SiO and HCO + show the same distribution 
of the high velocity gas. Here we only present the results from our HCO + imaging, which 
has better signal-to-noise and higher spatial resolution. 

HCO + shows a strong, compact, jet-like bipolar high velocity outflow at a pa of 147°± 5° 
(Fig 3). The symmetry axis of the outflow passes through the VLA source. The H 13 CN map 
shows that the surrounding cloud core is more dense to the north than to the south. This 
is also evident from the outflow, which is much more compact to the northwest than to the 
southeast. In the northwest the outflow is blue-shifted and terminates sharply at ~7" from 
the protostar. To the southeast the outflow is mostly red-shifted and well collimated (aspect 
ratio ~ 2:1) to ~16" from the protostar, where it expands to become a more wide angle flow, 
presumably breaking out of the cloud at ~ 30" from the protostar. The strongest blue-shifted 
emission is seen at the tip of the northwestern outflow lobe with outflow velocities extending 
to 27 km s _1 from the systemic cloud velocity. The southeastern outflow lobe shows both 
red- and blue-shifted emission with outflow velocities of ~24 km s -1 in the red and blue- 
shifted emission up to 18 km s~ x in the dense, well collimated part of the outflow. At the 
tip of the outflow, the outflow velocities are lower, presumably because the high velocity gas 
is too tenuous to be excited in HCO + . 

Since there is both blue- and red-shifted emission to the southeast , and little overlap 
between the outflow lobes at the star, the outflow must have fairly high inclination. In the 
following we assume an inclination of 60°, which is roughly consistent with the observed 
aspect ratio of the disk, i.e. we see the disk almost edge on. The northwestern blue-shifted 
outflow lobe is still confined by the surrounding cloud, while the southeast outflow appears 
to have broken through the dense cloud core. Since the extent of the outflow to the southeast 



- 6- 



is less well defined we derive outflow characteristics mainly from the northwest blue-shifted 
outflow lobe. 

The blue-shifted outflow lobe gives a dynamical time scale of ~2,000 yr, while the red- 
shifted emission gives ~10,000 yr, corrected for a 60° inclination. Note that the velocities 
are probably much higher in the red-shifted outflow, since it is expanding through a less 
dense cloud. We assume HCO + to be optically thin with an excitation temperature of 30 K, 
the same as Shepherd and Churchwell (1996) used in their study of outflows from high mass 
stars. If we further assume a normal HCO + abundance in the outflow (Section 4), we 
find a mass of 9 M Q for the northwest blue-shifted outflow and 8 M Q for the redshifted 
gas in the southeast outflow lobe. The mass weighted momentum flux, F = mv/td yn and 
the mechanical luminosity, L mec = mv 3 /2r, where r is the extent of the outflow, is 90 10~ 3 
M & kms~ l yr~ l , and 130 L Q for the blue outflow alone. The mass loss rate, based only on 
the confined blue-shifted outflow, could be as high as ~5 10 -3 M Q yr^ 1 . 

6. Discussion and Conclusions 

The new data we present in our letter reveal a well resolved elliptical source mapped in 
several optically thin lines as well as in 3.4 mm continuum. The peak of the line emission is 
centered within 1" - 2" of the VLA 6cm continuum, OH, and H 2 maser emission position. 
The 3.4 mm continuum emission has an extent similar to the elongated source seen in 
optically thin high-density tracers with a peak (hotspot) within 1" of the VLA and maser 
position. Follow-up observations at 1 mm in both line and continuum (Sandell, Wright, & 
Forster 2003, in preparation) show the disk more clearly and show that CH 3 CN J— 12— 
which traces hot gas, is centered within 0''2 of this position. FIR observations (Werner et 
al. 1979; Thronson and Harper 1979) show a FIR source with a luminosity of ~ 10 4 L , 
with a narrow FIR spectral energy distribution, strongly suggesting that the FIR luminosity 
is produced by a single central source. The FIR source has no near or mid-IR counterpart 
suggesting that it is heavily obscured. The FIR observations do not have enough spatial 
resolution or positional accuracy to confirm that the luminosity originates from the the 
OH/H2O masers, but all 1665 MHz OH masers require a luminous source for their excitation; 
therefore we conclude that the source powering the OH maser is the origin for the bulk of 
the observed FIR luminosity. 

The H 13 CN map (Fig 2) shows a clear velocity gradient along the major axis of the 
elliptical source. If we interpret this as a rotating disk gravitationally bound by the central 
source, we find that the center of symmetry is within < 2" of the maser source. Since this 
position also lies on the symmetry axis of the outflow, defined by our high spatial resolution 
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HCO + map, there is no doubt that this is the center of the mass and the protostar. From 
the observed velocity gradient we derive an enclosed mass of 40 M (Section 4). Since we 
have resolved the emission of the disk from the envelope, we can use the integrated emission 
of the disk to derive a mass. Assuming normal abundances and a kinetic temperature of 
35 K for the disk we find a mass ~ 400 M Q from the H 13 CN, and ~ 100 M Q from the 
3.4 mm continuum emission. Considering the uncertainties in molecular abundances and 
dust opacity, these estimates are in reasonable agreement. The mass of the central object 
and the disk are of the same magnitude. 

The bipolar HCO+ high velocity outflow (Section 5) may be driven by the rotating 
protostellar disk. The outflow is extremely young; the dynamical timescale from the blue- 
shifted outflow indicates an age of < 2000 yr, with an upper limit of 10,000 yr derived 
from the red-shifted outflow. Shepherd and Churchwell (1996) show that there is a clear 
correlation between mass loss rate and bolometric luminosity of the central source over a 
wide range of luminosities. If we use the correlation derived by Shepherd and Churchwell 
(1996), our observed mass loss rate corresponds to a luminosity of close to 10 5 L Q , while the 
observed luminosity is only 1O 4 L . Our observations therefore suggest an elevated mass loss 
rate in the very earliest stages of the formation of a massive star. Since we find the rotating 
accretion disk to be more massive than the central protostar, the disk will be highly unstable 
and is likely to feed both the outflow and the central accreting protostar. Model calculations 
by Yorke, Bodenheimer and Laughlin (1995) for a somewhat less massive central star (10 
M ), find that the accretion disk evolves into a thin flared disk with a ratio of M±/M disk ~ 
4, for all the models that they computed. Even though they ignored magnetic fields, these 
models should still give an idea of how the disk is expected to evolve. 

There is only one near-IR source in the vicinity of NGC7538S, IRS 11 (Fig 3). This 
IR source is associated with nebulosity and a knot of vibrationally excited H 2 emission, 
suggesting that it is most likely a young star. There are several additional H 2 knots in 
the vicinity of NGC7538S; two appear to be associated with the molecular outflow. The 
only other evidence for young objects in this region is an unexplained high velocity HCO + 
emission feature 10" southwest of the protostar, three vibrationally excited H 2 knots ~ 20" 
east of it, and an H 2 maser ~ 30" west of it. The high velocity HCO + could be part of 
a wide angle disk-wind or due to a young low luminosity source in the outskirt of the disk. 
This could explain the observed asymmetry in the continuum disk (more extended to the 
southwest), but it cannot produce the observed velocity gradient. The eastern H 2 knots and 
the western H 2 maser are unrelated to the protostellar source. 

It is therefore possible that the extended disk/envelope system harbours yet another 
protostar. Such a protostar would have much lower luminosity and mass, since the massive 
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protostar can explain most, if not all of the observed luminosity. Although high- mass stars 
usually form in clusters, all the evidence we have so far, suggest that NGC 7538 S is an 
isolated high mass star. If the central protostar is a single star, it may evolve into a late O 
or early B-star. 

We thank Dr. Osamu Kameya for sharing his unpublished VLA results with us. The 
BIMA array is operated by the Universities of California (Berkeley), Illinois, and Maryland 
with support from the National Science Foundation. 
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Fig. 1. — Color image of 3.4 mm continuum emission of NGC 7538 S overlaid with contours 
of H 13 CN J=l— emission integrated over the velocity range —70 to —40 km s _1 . The maps 
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are offset (-3"7,-0"5) from the 3.4 mm peak: a(2000.0) = 23 ft 13 m 44598, 5(2000.0) = + 
61° 26' 49''2. The intensity scale for the continuum is in Jy beam -1 . The contours for the 
H 13 CN emission go from —6 K km s _1 (dotted contours) and thereafter from 6 K km s _1 
with steps of 6 K km s _1 . The beam FWHM is plotted in blue. 
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Fig. 2. — Position velocity plot of the H 13 CN J=1^0 emission along the disk plane at pa 
= 60.0 °. Positive offsets are to the northeast. The velocity gradient is seen in all three 
hyperfine lines. The velocity scale is relative to the strongest hyperfine component, F = 2 - 
1, centered at a systemic velocity of ~ —56 km s _1 . The line drawing at the bottom of the 
plot shows the position of three hyperfine components and their relative intensity. The disk 
emission is red-shifted to the southwest and blue-shifted to the northeast (see text). 
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Fig. 3. — Contour map of blue- and red-shifted high velocity HCO + J = 1^0 integrated 
over 20 km s -1 overlaid on a greyscale image of integrated H 13 CN. The red circle marks the 
position of the protostar and the white line shows the symmetry axis of the outflow. The 
star symbol marks the position of IRS 11 and the filled squares are H 2 knots from Davis et 
al. (1998). The contours are logarithmic and start at 4.6 K km s" 1 for HCO + with a step of 
10°- 2N , where N= 1,2.. The beam FWHM is plotted in green. 



